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We introduce a positive Hermitian operator, the Fisher operator, and use it to examine a measurement process 
incorporating unitary dynamics and complete measurements. We develop the idea of information complement, 
the minimization of which establishes the optimal precision for a fixed input. The formalism demonstrates 
that, in general, the classical Fisher Information has the Hamiltonian semi-norm as an upper bound. This 
is achievable with a qubit probe and only projective measurements, and is independent of the true value of 
the estimated parameter In an interferometry context, we show that an optimal measurement scheme can be 
constructed from linear optics and photon counting, without recourse to generalised measurements or exotic 
unitaries outside of SU (2). 
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Scientists strive for an understanding of Nature by a physi- 
cal interaction introducing correlations between observer and 
observed, and the process is called measurement. Quantum 
mechanics exerts fundamental limitations on the precision of 
any measurement [l], yet for a quantum observable there is 
still a classical probability distribution over measurement out- 
comes. This distribution may depend on some real-valued 
system parameter 9 such as an interaction time or phase an- 
gle. It may have no associated Hermitian observable, nor be 
measurable directly - but its estimation may be the true goal 
of the measurement. Inferring 6 from frequencies of measure- 
ment outcomes is a conventional challenge in classical infor- 
mation theory with a well-established methodology |2|. From 
this classical underpinning we will show a quantum formalism 
emerges directly and without recourse to the quantum Fisher 
Information IJIIU or a linear error propagation model 13] 0|. 
Those approaches focus on the geometry of the state in Hilbert 
space and have a long distinguished history stretching at least 
as far back as 1976 |7|. In contrast, the formulation in this 
paper is self-contained and incorporates all three instrument 
components - input, dynamics and measurement choice - on 
an equal footing. Our conclusions rely on no other results or 
assumptions than those presented herein. 

Consider a maximal test 1 8 1 having outcomes labelled k and 
an associated probability distribution P{0) = {pa;(^)} that 
depends on a continuous real parameter 9. Classical Fisher 
information |9| is defined as 

:7(0)=^Pfe((?) (^J^lnpfc(0)^ . (1) 

This functional is a measure of the information contained in 
the distribution P{9) about the parameter 9 It provides 
the sole distance metric in probability space ifTOl and the scal- 
ing factor in local distinguishability ifTTl . An explicit lower 
bound on the standard error of an unbiased estimate 9 on the 
true phase 9 is given by the reciprocal of the Fisher infor- 
mation, (39)'^ > 1/Ji9), called the Cramer-Rao bound |9|. 
Therefore, for optimal precision one maximizes the Fisher in- 
formation 1121. 



Let's take the case of the measurement observable M and 
outcomes {rrik} associated with distribution {pk} applied to 
a quantum system that previously evolved from a known ini- 
tial state \tpo) under the dynamics of some Hamiltonian H 
for time 9. The Schrodinger equation governs the dynamics 
idgltpe) = H\tpe), (where h = 1) and the time evolution is 
expHcidy \^g) = exp{—iH9}\^o) . 

The parameter estimation task involves an inference of 
the time-like variable 9 from the measurement distribution 
{pk{9)}. Writing the spectral decomposition of the maxi- 
mal measurement as M = J2k''^k\k){k\ then complex am- 
plitudes {k\ipg) = rkexp{i(j)k} lead to probabilities pk = 
{k\il;g){tpg\k) = r| , where it is understood that {pk,rk} G 
[0, 1] and (j)k € [0, 27r) are all real-valued functions of param- 
eter 9. Replacing pk with in Eq.Q gives: 

J{9)^iJ2^l, (2) 

fc 

which we can now use to find an operator expression for the 
Fisher information. Differentiating rf. gives 

2rkrk = {k\^e){Mk) + {k\iPg) {iPe\k) . (3) 

Now, {k\tpg) = dg{rke^'>"') = e"*"'{f + ir^k) = 
I (fcl-^e) |e*('^'=+'^'=^ where we define the inclination of the ve- 
locity vector: 

Tfe = arg(fc|i/'e) - arg(fc|'09) = tan"^(rfc0fe/ffe) . (4) 

Eq. ([3| yields ffe = cos Tk\{k\%l)g)\. From the Schrodinger 
equation we also have = {k\H\'il)g). Substituting this 

expression and squaring gives 

tI = cos" Tk{MH\k){k\H\i:g). (5) 

Summing over all the measurement outcomes 'fc' 
gives the Fisher information using Eq.(|2|: J{9) = 
4 X]fe '^os'^ Tfe(^e|_ff |fc)(fc|_ff l?/'^). Therefore, we can de- 
fine a non-linear, positive and hermitian Fisher operator Fg 



2 



that is diagonal in the measurement basis: 

Fe = 4^cos2 rfe|fc)(fc| = 4^ cv,,fc|fc)(fc| , (6) 

k k 

such that the Fisher information is then 



J{d) = {MHFgH\i^e) = Tr \^e){M HFgH 



(7) 



This can be re-expressed as an expectation value with re- 
spect to the probe state, the input |V'o)^ by unitarily trans- 
forming the Fisher operator: J^{0) = {ip()\H^gH\tljo) = 
\ijQ){iJo\H^eH , where 

" Fee-'"' = AY,C4,,k\k'){k'\ (8) 



Tr 



= e 



.H0 



is the unitarily transformed Fisher operator, cf Eq. (|6|. Note 
that, due to the non-linear nature of Fg, a basis transformation 
|fc) I— > |fc') would give a different result: 



k 



(9) 



not equivalent to ^g, since c^,k' 7^ c^,k generally. 

Fixed Probe Optimization: The completeness of the 
measurement basis provides a resolution of the identity 
\k) {k\ ~ ^iid therefore Eq.(|7]) becomes: 



-^sin^Tfc {i:g\H\k){k\H\i:g) 



^{rk(l>k) 



(10) 



where the form in the second line was previously presented 
in ifTTl . With reference to Eq.(|2| it is interesting to note 
that (-ff^) is a sum of translational and rotational kinetic en- 
ergy terms. Here it may be taken with respect to IV'o)- since 
[H'^,exp{iH9}] = 0. In Eq.([TO| the subtracted term, we will 
call it the 'information complement' JC, is positive and can 
only reduce J {6). We look for a measurement basis {|fc)} 
and associated set {rk,(f)k} that minimizes this information 
complement for a fixed input It/jq)- Finding derivatives of (f)k 
by writing 54>k = arg(/c|V'e+5e) - a.rg{k\ipg) produces 



rk(Pk 



^ \Hkj\riC0s{4>i - 4>k + ^k,i) = Ak , (11) 



where {k\H\l) = Hk^i = \Hk,i\expi^k,i, and therefore 

/C-^(rfe0fc)2=^^^ (12) 

k k 

By comparison, expanding {H) in the \k) basis gives 
{IF) = ^rirfc|i/fc,i|rzcos(0; - (j)k + ^k,i) = ^ »'fc A 

k,l k 

(13) 




FIG. 1: Left: The optimal measurement scheme for any Hamiltonian 
a may be restricted to the qubit subspace of its extremal eigenvec- 
tors |A|) and On the Bloch sphere: the diagonal basis of fl 
defines the z-axis, and the optimal measurement projectors \ k±) the 
a;-axis. An optimal probe state IV'opt) lies anywhere in the equato- 
rial plane. Right: The vector construction for finding c^.^ in Eq.l[22} 
with A> 1. 



At the minimum 



— - V— (A) - 



(14) 



and this calculation leads to 

^ ^ _^ Efc A|gfc,p| sin(0p -(t>k + ^k..p) Q 

^ ^ J2k'^k\Hk.p\sm{cj)p- (l)k + ilk,p) ^ ^' 

(15) 

at the minimum of /C. Remember, this is equivalent to max- 
imizing J' over {\k)} for a fixed IV'o)- At this minimum we 
can prove the relation that 

IC - (H)^ 

=Y.^l-{Y.-kAk)' 

k k 

= Y.'lBl^{Y,rlBkf 

k k 

= A2s>0, (16) 



since {r\} = {pfe} is a probability distribution. Comparing 
again with Eq.([TO]l it follows directly that 

J(6i)/4<A2iJ (17) 

for a fixed input ItAo)- Now we will show that this bound is 
saturated by a particular qubit input state. 

Optimizing for a Qubit: For probe states that are eigenstate 
of the Hamiltonian then all rk ^ Q and J ~ Q from Eq.(|2]) 
because any eigenstate |A) of _ff only gains a phase during 
its evolution. Thus rk{9) |(fc|e^-^^|A)| = Kfele^^^lA)] = 
|(fc|A)| = rfc(O), and fk = 0. For optimality over all j-^o) and 
{\k)} the input state must thus be a superposition of at least 
two Hamiltonian eigenvectors. 



I^/iq) ^ COS + e'^ sin 7] A2) 



(18) 



Starting with a qubit probe state of the form of Eq. ( [T8] l and 
some measurement basis pair {[/ci), 1^2)} spanning the same 
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C2 as{|Ai),|A2)}: 



with A= A2/A1 and 



\ki) = cosa|Ai) + sina|A2) 
1/02) — — sinalAi) + cosa|A2) 



(19) 



Here it has been chosen that the measurement basis defines 
the X axis on the Bloch sphere of FIG[T] hence the real coef- 
ficients for |fci,2)- By confining the measurement basis vec- 
tors {|fci,), 1^2)} to the same space as {|Ai), IA2)} then one 
can restrict interest to the component of <tg within this qubit 
space: 



$e = e'^^{c^,fci|fci)(/ci| +c^,fc2|A:2)(fc2| + ■ ■ -je 



~iHe 



-j(A2-Ai)fl 



c-4Kki I +i(\2-\i)e 



_g-«(A2-Ai)e)g^ 



+ C^,k2 [ _g+j(A2-Ai)e)gg ^2 



where c (s) is cos a (sin a). We ignore elements of $£/ 
that project onto the remaining Hilbert space, orthogonal to 
{|Ai), IA2)}. The coefficients cos^ have a geometric inter- 
pretation by mapping (D i-^ R^: 



Tfc = cos 



1 (fc|^e).(fc|^g) 

m^e)\\{km\ 

I {v + Aw).(v + w) 
\v + Aw\ \v + w\ 



{k\H\4,e).{k\^e) 
mH\^e)\mrPe)\ 

(20) 



ui = cos(a) cos(7), zwi = sin(Q!) sin(7)e**''^ 

V2 = — sin(a) cos(7), W2 = cos(a) sin(7)e*^'^~*-'^^^'^^^^' . 

The components in the measurement basis are: 

(fci,2|-ffe-^^Vo) = e-*^^^Ai(t/i^2 + Aw^^2) (21) 

and the included angle between vectors {v + w) and {v + Aw) 
is Tk — 7r/2 + arccos{Ai/|Ai|}, because Ai may be negative. 
Now we have explicit expressions for 

c^,fe = cos^ Tfe = (1 - sin^ Tfe) = 1- 

{AR^ + 1) cos(/3)i?+ 1)^ 

(i?2+2(-l)fe-i cos(/3)i?+l)(^2^2 4_2(-l)fe-iAcos(/3)i?+l) 

(22) 



with 

{i?l,i?2}^ 



/Kl 


KM 




I l^^il 


'l«2| J 





tan(7) \ 
tan (a) J 



(23) 

and (3 — arccos{wl.wi/|wi||wi|} = X ^ (-^2 — Ai)6'. For 
angles {a, /?, 7} expectation value (■(/'o|-H^*&9i?|-0o) gives 



J(a,/3,7) = 



-4(Ai-A2)2s2[2a]s2[27]s' 



(c[2(a - 7)] + c[2(a + 7)] + 2c [/?] s[2a]s[27] - 2)(c[2(a - 7)] + c[2{a + 7)] + 2c [[3] s[2a]s[27] + 2) 



(24) 



writing sin as 's' and cos as 'c'. This function is optimized by 
angles {a, 7} 1-^ 7r/4, independent of the value of /? and giv- 
ing a saturable bound: (V'ol-ff'&e-ff iV'o) < (Ai — A2)^. This 
is the upper bound on the Fisher information for any superpo- 
sition of two eigenstates of the Hamiltonian. It is saturated by 
a probe state 



|^op,) = (|Ai) + e^^|A2))/V2 



(25) 



where x € [0, 27r), see FIG{T| The result a 1-^ 7r/4 dictates 
an optimal measurement scheme with components: 



(|Ai)±|A2))/V2 



(26) 



also depicted in FIG{T| All other basis elements G {\k)} span 
an orthogonal subspace. 

Generalization to Higher Dimensions: It is clear that for a 
given Hamiltonian H, the maximal Fisher information is then 
bounded from below by (A| — X^f = where A| (A^) 

is the max (min) eigenvalue of H, and \ \H\\ — (A| — A|) is 



the operator seminorm of the Hamiltonian ||4| . A key property 
of the seminorm is that it gives an achievable upper bound to 
the variance: 



\H\V > AA^H. 



(27) 



This allows a connection to be made between the qubit re- 
sult with that for a fixed |V^o) - we saw earher for a fixed 
input that J' < iA^H. Therefore the qubit maximum vari- 
ance state must be the universally optimal state over the full 
Hilbert space, as it saturates its variance bound = 
The general result is 



\i^o)A\k)} 



J 



\H\ 



(28) 



A corollary of Eq.(|28)l is that no greater number of superposed 
energy eigenstates can improve on the Fisher information pro- 
vided by the maximum variance state (|A|) + e^^\Xi)) /y/2. 

Photons: In a Mach-Zehnder interferometer (MZ) the esti- 
mated parameter 6 is simply the phase difference between the 
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two interferometer paths. Spin jy plays the role of Hamilto- 
nian, and the Casimir operator ,P is equivalently n/2(n/2+l) 
where h is the total photon number operator (So j = n/2 for 
states of fixed n.) Eigen-equations are J^|j, to)^ — j{j + 
l)\j,m)i and Ji\j,m).i = m\j,m)i fori € {x,y,z}. The 
maximum variance state for the MZ is 

IV'opt) = +J% + e'^^b, , (29) 

i.e. a NOON state HI] [HI, rotated by 7r/2 around the x 
axis, with arbitray phase x- Eq.(|26]) indicates that the opti- 
mal measurement basis has two elements |fc±) = (Ij, +j)y ± 
exp(i^) |j, —j)y)/ v^- Interestingly, the optimal measurement 
basis is not unique (in this case) and maximal precision is also 
recovered by {|/c)} corresponding to the eigenbasis of J 2, a 
measurement of photon number difference between the two 
interferometer modes ifTTl . Therefore, for a lossless MZ the 
precision limit is saturated only by NOON probes of Eq.(|29|, 
and this is achievable with simple projection measurements, 
i.e. without generalised measurements - compare lfT4l [TSl . 
These optimal projection measurements may be performed 
by linear optics and photon counting alone. (For any pho- 
ton number space of n photons, of the full group of Unitary 
operations, SU(n + 1) , only the SU(2) subset is needed.) 
The Cramer-Rao inequality in the context of interferometry 
thus gives a precision known as the 'Heisenberg' limit, iTTSl : 
(-^^h)' > 1/X,ax = 1/(At - h? = l/4j2 = l/n\ 
Phase States As an example of a fixed probe, take the phase 

state mini 

(parameterized by C G R)- In a MZ the evolution is a transla- 
tion of the parameter, e^*^"^" |j, C,) — |j, ( — 6). These states 
are of interest because it is expected that their precision scales 
J' cx j^, close to the Heisenberg limit. They have the peculiar 
feature that, for any m, 

arg ,{j,m\d\j,e) = -j7r/2 , V 6 G 50(2j + 1) , (31) 

due to the properties fTSl of the Wigner rotation elements 
2O, Wile"*"''' ^ I j, 77^2)2 G ]R. It follows that choosing \k) = 
0\j, m)z makes 0^ independent of 9, i.e. = i/ife = /C = 0. 
Therefore Eq.([TO| gives J^,^ ^ 4(j2) = 4j(j + l)/3, 
with the desired scaling. 

Summary and Outlook: The formalism developed here in- 
corporates all aspects of quantum parameter estimation ex- 
plicitly; probe IV-'o). dynamics H, and measurement, 
clarifying how precision is determined by the interplay of all 
three. We introduced the information complement JC, the min- 
imization of which allowed the optimal measurement to be 
found for a fixed probe. The greatest precision was found in 
the qubit subspace of the maximal variance input state. This 
ultimate precision is completely defined by the difference of 



the extremal energy eigenvalues - no deeper dynamical struc- 
ture is relevant, nor is the dimension of the Hilbert space. 

Despite its many advantages (such as its easy extension 
to mixed states), if precision limits had been found using 
the quantum Fisher Information |3j i4j|, nothing could be 
learned about the relative performance of various measure- 
ment schemes, as that approach assumes an unknown optimal 
measurement. But for specific dynamical processes (e.g. pho- 
ton interferometry) there exist real physical restrictions on the 
types of probe and measurement that may be employed. For a 
restricted set that excludes an optimal measurement, maximiz- 
ing the classical Fisher information over the available mea- 
surements will allow the locally optimal probe to be recov- 
ered. 

In future, it may prove fruitful to develop the Fisher Oper- 
ator approach to incorporate evolution of mixed states, gov- 
erned by completely positive maps and generalised measure- 
ments. 
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